Supersolid phases of a doped valence-bond quantum antiferromagnet: 
Evidence for a coexisting superconducting order parameter 
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Motivated by numerical evidence of the valence bond groundstate of the two-dimensional Heisen- 
berg pyrochlore lattice, we argue using a t-J model that it evolves under doping into novel phases 
characterized by superconductivity coexisting with the underlying valence-bond solid order. A 
fermionic mean-field theory supplemented by exact diagonalization results provide strong arguments 
in favor of the stability of such supersolid phases. The resemblance with modulated superconducting 
patterns in high-T c cuprates as well as possible relevance to frustrated noncuprate superconductors 
such as spinels and pyrochlores is discussed. 
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The so-called supersolid (SS) [l|, an exotic state of 
matter that breaks both the spatial lattice translation 
symmetry and the internal E/(l) symmetry, associated 
with particle number conservation, has attracted much 
attention following its possible discovery in the solid 
phase of 4 He [2]. On the one hand, studies of various 
models with hard-core bosons on a frustrated triangular 
lattice demonstrate that bosons may condense forming a 
superfluid on top of the crystalline background 3 . More- 
over, there has been a convergence of agreement concern- 
ing nonuniform condensation of magnons in spin models 
with external magnetic field suggesting an exciting 
possibility of the formation of a spin SS in real quantum 
magnets. On the other hand, high-T c superconductiv- 
ity potentially described by the resonating valence bond 
(RVB) state on the square lattice with a superconduct- 
ing (SC) parameter of d- wave orbital symmetry turns out 
to be robust against spontaneous breaking of translation 
symmetry [El, Q. Hence, charge order coexisting with 
modulated SC order in Ca2-a:Naa;Cu02Cl2 revealed by 
recent scanning tunneling microscopy [3] may well have 
an extrinsic origin ||. 

It is therefore legitimate to question the relevance of a 
fermionic analog of the bosonic SS and its possible stabi- 
lization by geometrical magnetic frustration. In this re- 
spect, among of two-dimensional structures, the quantum 
antiferromagnet on the checkerboard lattice (see Fig. [T]) 
provides a priori conducive conditions to the appearance 
of a fermionic SS upon doping. Let us briefly summa- 
rize remarkable properties of the checkerboard lattice: 
(i) Made out of corner-sharing plaquettes with crossed 
bonds, it might be considered as a projection of a more re- 
alistic corner-sharing tetrahedra pyrochlore lattice known 
to host superconductivity, e.g., spinel LiTi2 04, as well 
as a- Cd2Re207 and /3-pyrochlore KOS2O6 @; (ii) Ex- 
act diagonalization (ED) [lOj and other intensive studies 
of the spin- 1/2 Heisenberg model on the checkerboard 
lattice [ll| suggest that its ground state is twofold de- 
generate being a product of plaquette singlets. It breaks 



translation symmetry while preserving the C& v rotation 
symmetry around the center of the void plaquette and 
as such it is a good example of a valence bond crystal 
(VBC) [12|; (iii) Lastly, particularly important sugges- 
tion for supersolidity in this case comes from a recent 
discovery of hole pairing [l3| , since alone the existence of 
a spin gap at half-filling is not sufficient to guaranty SC 
behavior at finite doping 
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In this Letter we show that a fermionic valence bond 
supersolid (VBSS), characterized by the coexistence of 
the SC order parameter and solid component of spin sin- 
glet order, can appear by doping the VBC characteriz- 
ing the Mott phase of some frustrated lattices like the 
checkerboard lattice. Interestingly, such a scenario has 
close correspondence with the RVB theory in which the 
translationally invariant <i-wave superconductor emerges 
out of the spin-liquid Mott phase. The alternative that 
the VBC is destroyed upon doping leading to an exotic 
d+id and d+is superconductor [15| with broken time- 





FIG. 1: (color online) Schematic pattern of dr^ i7r )-VBSS (a) 
and Sjv^j-VBSS (b) phases. In the checkerboard lattice, di- 
agonal bonds are on every second plaquettes. Rotation sym- 
metry of the pairing amplitude Ay reduces the number of 
inequivalent bonds down to three. The line widths of these 
bonds labelled by a ( = "if running from 1 to 3) are propor- 
tional to the magnitude of the corresponding spin-spin corre- 
lations (for t! It = —1). Different types of those lines (solid 
and dashed) correspond to opposite signs of A a . 
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reversal symmetry is being reconsidered [16J. In con- 
trast, we find a very rich phase diagram with various 
VBSS phases as summarized in Fig. O These phases dif- 
fer from their bosonic analog not only in the way lattice 
symmetry is broken, exhibiting bond modulation rather 
than charge order but, in addition, from the fact that 
off-diagonal and solid orders involve different (charge and 
spin) degrees of freedom. 

To address this intriguing issue, we consider a t-J 
model on the checkerboard lattice, at doping x = 1 — n, 
n being the electron density (per site): 



s-RVB 



n = - 



E 



(i) 



(ij) 



where the sums run over all bonds of the checkerboard 
lattice. For sake of generality, we allow for different hop- 
ping amplitudes and antiferromagnetic superexchangc 
couplings on the diagonal bonds, — t' and J,j = J\ 
and on the vertical and horizontal bonds, Uj — t and 
Jij = J. Note that under electron-hole transforma- 
tion t' changes into —t' so that hole and electron doping 
correspond to opposite signs of t' (while t > can al- 
ways be assumed). Hereafter, we assume a typical value 
t/J = 3 and set J'/ J = (t'/t) 2 which follows from the 
superexchange relation J — 4t 2 /U valid in the large U 
limit of the Hubbard model. Next, we replace local con- 
straints that restrict electron creation operators cj CT to 
the subspace with no doubly occupied sites by statisti- 
cal Gutzwiller weights, while decoupling in both particle- 
hole and particle-particle channels yields a renormalized 
mean- field theory (RMFT) Hamiltonian [l7j], 



H = 



~ E 9ijtij(4,a C j,<T + h ' c -) -ME™*' 



7 E 9ijJij(Xji4,a c j,<r + h-c.-Xij) 



- \ E 4 J v( A AA-« + h - c - - i a ^i 2 )' ( 2 ) 

with the Bogoliubov-de Gennes self-consistency condi- 
tions for the bond- Xji = ( c ] a c i.u) (assumed here to 
be real) and pair-order Aji = k r c^J = (c i _ a Cj a ) 
parameters in the unprojected state that we solve on 
a 48x48 cluster at low temperature (3 J — 200. Note 
that since we restrict our study to strong frustration 
(\t'\ ~ t), we do not consider here long-range antiferro- 
magnetic order, known to coexist with superconductivity 
at low doping on the square lattice l£j . Finally, to allow 
for small nonuniform charge modulations (if any), the 
Gutzwiller weights have been expressed in terms of local 
doped hole densities Xj = 1 — ^2 a (c\a c ia) as 9ij = \J z i z j 
and gfj = (2 — Zi)(2 — Zj) with Zi — 2xi/ (I + Xi). 

We mostly concentrate here on the t' < case with 
emergent supersolidity. Guided by the earlier studies at 




FIG. 2: (color online) RMFT phase diagram of the generalized 
t-J model on the checkerboard lattice with t' > (left) t' < 
(right) while J'/J = {t'/t) 2 . Solid (dashed) lines show first- 
(second-order) phase transition, respectively; striped area cor- 
responds to phase separation. 



half-filling [H1E3, 

in addition to usual d- and s-RVB su- 
perconductivity we also allow for coexisting VBC order in 
the mean-field equations. The d-wave symmetry 17| is a 
natural choice compatible with the required C^ v rotation 
symmetry of |A,j| around the centers of void plaquettes 
(bonds labeled by I in Fig. [T]). It yields at half-filling 
a strong insulating VBC with decoupled void plaquette 
singlets (i.e. \a = A a = on the a = 2, 3 bonds) which 
provides the lowest magnetic energy Ej/J = —3/4 in 
good agreement with numerics (lfj| . Upon finite doping 
x, an additional A 3 component arises along the diagonals 
as depicted in Fig.[]7a). Unit cell doubling together with 
d-wave symmetry implies a (tt, tt) phase modulation of 
A3 so that this phase is referred to as d^^-VBSS. In 
fact, away from half-filling, another SS state, labeled by 
S( W)7r )-VBSS and shown in Fig. [DJb), is competing with 
the d( ff-7r )-VBSS. In this phase, all neighboring void pla- 
quettes have opposite signs of the pairing parameters and 
s-wave rotation symmetry on every void plaquette. 

The spontaneous lattice symmetry breaking of the dis- 
covered VBSS phases is seen in Fig.[3Ja-d) showing bond 
charge hopping T a — 2g t a x a an d spin-spin correlations 
S a = — 1<?„ (xi + |AJ 2 ) (corresponding to the two terms 
of Eq. (J]), respectively) on the three inequivalent bonds 
defined in Fig. [TJ The s^ ^-VBSS is characterized by a 
weaker modulation of S a as compared to its d-wave coun- 
terpart. Even though the d-wave RVB order is known 
to yield a better kinetic energy than the s-wave RVB 
one 11711 ■ the s^^-VBSS has better total kinetic energy 
E t [16( due to the absence of phase modulation in the 
pairing A 3 . Next, in Fig.[S[e,f), we present the SC order 
parameter components A^ c which are related in RMFT 
to pairing amplitudes by A^ c = <7oA a and hence van- 



ish as expected at half-filling [I7J. Remarkably, all of 



them show a non-monotonic dome shape as a function of 
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FIG. 3: (color online) (a,b) Spin-spin correlations S a and (c,d) 
bond charge hoppings T a on the three inequivalent bonds 
of the d( XjW )-VBSS (left) and s^^j-VBSS (right) phase for 
t' jt = — 1. Panels (e,f) show the SC order parameters A^ 
found to be enhanced with respect to Aq C = A?£ = Af c on 
a simple square lattice [if = 0), see also Ref. hjl. Shaded 
areas denote the region(s) of relative stability of each phase. 



doping. Although the Gutzwiller RMFT is a zero tem- 
perature scheme, we expect that SC properties will be 
boosted to higher temperatures compared to the d-RVB 
phase with t' = 0. Note that VBC and SC orders can 
give rise to separate critical T c . Finally, the low-energy 
fermionic quasiparticlc band structures of the SS phases 
are shown in Fig.SJ Interestingly, in the d( 7r 7r )-VBSS, the 
A3 s-wave component on the diagonal bonds suppresses 
the usual zero-energy Dirac cone of the d-wave supercon- 
ductor along its (0, 0)-(tt/2, tt/2) nodal direction. Upon 
further doping, a gapless spectrum is recovered in the d- 
RVB phase with A 3 = due to melting of the valence 
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FIG. 4: (color online) Quasiparticle band structure found for 
t'/t = -1 at x = 1/4: (a) d^^-VBSS and (b) s^^-VBSS. 
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FIG. 5: (color online) Maxwell construction (see text) for 
t' < (a) and t' > (b). From top to bottom: \t'/t\=0.7, 
0.8, 0.9, 1.0, 1.05, and 1.1. The arrows indicate small kinks 
due to the first order transition expected if PS (cf. Fig. [2]| is 
not considered. 



bond order. In contrast, the energy gap of the s^.tt)- 
VBSS is less modulated as expected for isotropic pairing. 

It is known that doping a Mott insulator can lead to 
phase separation (PS), i.e., a macroscopic segregation 
into two phases, commonly the undoped phase and a 
phase with a finite doping 19]. Hence, we have per- 
formed the standard Maxwell construction and calcu- 
lated the slope of the ground state energy density S (x) = 
[E(x)—E(0)]/x as a function of doping x. Technically, PS 
is signaled by a negative curvature of E(x) when x — > 
(i.e. a negative compressibility) and a minimum of S(x) 
at a given doping x c (cf. Fig. [5|). We expect that be- 
low x c the system phase separates into a two-component 
mixture of the undoped (x=0) VBC and a doped phase 
(see the phase diagram) with x — x c . Note that other 
phase separated mixtures which could in principle ap- 
pear in the vicinity of some first-order transitions within 
the < x < x c region (seen e.g. as kinks in Fig. 5) are 
here prevented by the PS between the above-mentioned 
"extremal" phases. 

We now briefly discuss the case t' > 0. First, un- 
der doping, the emergent d( 7r 7r )-VBSS is now much more 
fragile w.r.t. the melting of the VBC order. Secondly, a 
smaller range of stability is found for the d-RVB phase 
compared to the previous t' < case consistently with 
the strong electron-hole asymmetry of pairing found in 
numerics [11 El . On top, as deduced from the Maxwell 
construction shown in Fig. [5jb), these narrow VBSS and 
c?-RVB regions appear to be unstable towards PS between 
a hole- free VBC and a homogeneous metallic Fermi liquid 
(FL) with x = x c ~ 0.15 at t' = t. 

Our findings are summarized in a global RMFT phase 
diagram in the (x, \t'/t\) plane reported in Fig. [2j It em- 
phasizes the crucial role of the frustration in stabilizing 
the VBSS states for t' < 0. Indeed, increasing frustra- 
tion strength \t'/t\ shifts the second-order phase tran- 
sition to the nearby d-wave RVB phase towards higher 
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doping level. A region of S( 7rj7r )-VBSS is squeezed into a 
narrow range 0.94 < \t' Jt\ < 1.07 limited below (above) 
by a first-order phase transition line to the ci^^-VBSS 
(s-RVB) phase. In fact, for maximum frustration, both 
d/ ni „y and s^^-VBSS states extend up to quarter-filling 
x = 1/2 where they merge to a VBC similar to the one 
at half- filling (albeit with no pairing, i.e., A a = 0). Note 
that in the whole PS region, the role of the so far ne- 
glected long-range Coulomb repulsion will become im- 
portant and lead to microscopic (possibly stripe-like) PS 
instead of true macroscopic PS. Also, although quantum 
fluctuations can play an important role, here we can nev- 
ertheless provide simple additional arguments supporting 
our mean-field phase diagram. First of all, the large spin 
gap (estimated by ED of a 32-site cluster to be around 
0.69 J [l(|) is expected to survive at least in the vicin- 
ity of half- filling. Next, we have found that the two-hole 
binding energy estimated by ED of the same cluster is 
substantial even in the J/t — > limit where it scales as 
—4.4396 J (with |£'|=t), confirming strong pairing. Fi- 
nally, we have checked that below \t'/t\ = 0.4, the VBSS 
melts towards phase-separated d-RVB state in agreement 
with Variational Monte Carlo simulations for t' = 21 1 . 

Experimentally, examples where spatial frustration 
and superconductivity coexist are KOS2O6 and LiTi2 04. 
Although electron correlations are usually associated to 
nodal pairing, recent experiments show the evidence for 
fully gapped superconductivity in both compounds [2^ |. 
so conventional phonon-mediated pairing has been advo- 
cated. However, our results indicate that, on a strongly 
frustrated lattice, electron correlations can also lead to 
gapful pairing and a modulated structure. We also note 
that, despite interesting similarities, the C^-symmetric 
VBSS states differ from the unidirectional C^v SS found 
on the square lattice of Ca2- a; Na : i;Cu02Cl2 [7j by the 
absence of spatial modulation in the local hole density. 

In summary, we have shown that the interplay between 
electron correlations and geometrical frustration can sta- 
bilize novel states of matter exhibiting microscopic coex- 
istence of superconductivity and spin dimer-crystalline 
order. In contrast to the inhomogeneous superconductor 
of the frustrated Shastry-Sutherland lattice which is dis- 
connected from the half-filled dimer ground state [23j |. 
fermionic supersolids emerge naturally upon doping a 
parent Mott insulator with spontaneous spin dimer or- 
der. We have identified on the checkerboard lattice two 
different competing supersolids with either s-wave or d- 
wave pairing symmetry w.r.t. the centers of the void pla- 
quettes. We believe that our findings open a new route 
to search for new exotic frustrated superconductors. 
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Supplementary data to the article: Supersolid 
phases of a doped valence-bond quantum 
antiferromagnet: Evidence for coexisting 
superconducting order 

In our Letter, we have provided strong arguments in 
favor of the formation of a fermionic valence bond su- 
persolid (VBSS) on the checkerboard lattice. It is a 
new alternative scenario to the one in which the va- 
lence bond crystal (VBC) known to set in at half-filling 
ToL 11 1 is destroyed upon doping while frustration pro- 
duces an exotic d+id and d+is superconductivity with 
broken time-reversal symmetry [151 ] . Hence, we have 
applied the renormalized mean-field theory (RMFT) to 
study the stability of various VBSS structures on equal 
footing with the d+id and d+is superconductors that 
preserve the translation symmetry of the lattice. The 
results for two representative cases: half-fillling x = 
and x = 1/4 doping (the latter chosen to be outside the 
expected phase separated region) are summarized in Ta- 
ble HI They clearly show that the half-filled Mott insu- 
lating ground state corresponds to the VBC with d-wave 
(spin) pairing amplitude and has substantially better en- 
ergy as compared to the next low-energy d+id order. Un- 
der doping it evolves into d( W)7r \-VBSS (in fact, thermo- 
dynamically unstable towards phase separation below a 

1 Similar "plaquette" SC state appears on the Shastry-Sutherland 
lattice, see B.-J. Yang, Y. B. Kim, J. Yu, and K. Park, Phys. 



small critical doping). Further doping results in a first- 
order phase transition towards a new s-wave supercon- 
ducting ground state which is found to coexist with (spin 
dimer) crystalline order forming thus the S( Wj7r )-VBSS. 
Note also that around x = 1/4 a simple homogeneous 
s (7r,7r)"RVB phase 1 with the same pairing symmetry as 
the s^^-VBSS (but with |Ai| = |A 2 | and A 3 = 0) also 
becomes very competitive. 



TABLE I: Site-normalized RMFT energy E and kinetic en- 
ergy Et of various wave functions considered in our Letter ob- 
tained for t' /t — — 1. Pure s-RVB and d+is phases converge 
at half-filling to a homogeneous Fermi liquid; d+id phase is 
unstable at x = 1/4 towards d-RVB. 





x = 


x = 


1/4 


phase 


E/J 


E/J 


E t /J 


s-RVB 




-2.0848 


-1.7928 


d+is 




-2.0868 


-1.7859 


S(^,^)-RVB 


-0.5944 


-2.1044 


-1.7964 


s (7r>7r) -VBSS/VBC 


-0.6134 


-2.1062 


-1.7945 


d-RVB 


-0.6884 


-2.0776 


-1.7752 


d+id 


-0.6886 






d (7ri7r) -VBSS/VBC 


-0.7500 


-2.0980 


-1.7729 
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